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Cation intermixing at the n-type polar LaGa0 3 /SrTi0 3 (001) interface is investigated by first principles 
calculations. Ti<s>Ga, SroLa, and SrTi<*LaGa intermixing are studied in comparison to each other, with a 
focus on the interface stability. We demonstrate in which cases intermixing is energetically favorable as 
compared to a clean interface. A depopulation of the Ti "id^ orbitals under cation intermixing is found, 
reflecting a complete suppression of the two-dimensional electron gas present at the clean interface. 

The high mobility (~ 10 3 cm'W 1 at 4.2 K) two-dimensional electron gas at the (LaO) + /(Ti0 2 )° n-type 
interface 1 between the band insulators LaA10 3 and SrTi0 3 (STO) is in the focus of interest since quite some 
time 2 4 . It has been proposed that several mechanisms such as the polar discontinuity 5 , structural relaxation 6 , 
interdiffusion of cations 711 , and oxygen vacancies 1214 , have to be taken into account to explain the interface 
conductivity. Of fundamental interest is also the thermodynamic stability of the interface due to the possibility of 
atomic intermixing in the interface region. Experimentally, cation intermixing of La and Sr at the LaA10 3 /STO 
interface has been observed, accompanied by considerable lattice deformations 7 . It is predicted that the interface 
polarization can be stabilized by the atomic intermixing 1517 . 

Recently, a detail analysis of different configurations of cation intermixing at the n-type LaA10 3 /STO interface 
has been performed 18 . It has been observed that the energetically most stable configuration is obtained when the 
interface dipole is cancelled out as a consequence of the cation intermixing. Little further attention has been paid 
to cation intermixing by first principles theory because of the heavy calculational efforts coming along with the 
large supercells required for such investigations. However, it has been found that intersite cation La/Sr disorder at 
the n-type LaTi0 3 /STO interface is energetically favorable as compared to an ideal interface, while the electronic 
properties of the system remain almost the same 19 . On the other hand, the metallic states at the DySc0 3 /STO 
interface 20 are suppressed for Dy<t>Sr and ScoTi intermixing 21 . It has been observed that the dipole energy at the 
(LaO) + /(Ti0 2 ) 0 interface increases when charge is transferred from the LaO layer to the adjacent STO unit cells 5 . 
This energy can be reduced by intermixing of the Sr and La cations to produce an electric field which compensates 
the interface dipole. 

In this work, we study TioGa, SroLa, and SrTi<s>LaGa cation intermixing at the n-type (LaO) l 7(Ti0 2 ) 0 
interface between LaGa0 3 (LGO) and STO. Formation of a two-dimensional electron gas at this interface has 
been reported experimentally and the effect of O vacancies on the conductivity has been addressed in Refs. 22-24. 
Metallic states appear for both the n andp-type LGO/STO interfaces, where the metalliciry of the n-type interface 
is enhanced by O vacancies. On the other hand, the p-type interface turns gradually from a hole doped into an 
electron doped state for strong O deficiency 25 . Polar distortion effects have been studied theoretically in Ref. 26. 



Results 

The bulk electronic structures of LGO and STO are well understood. Since the lattice mismatch between the 
component materials is small, the structural relaxation at the interface does not play a dominating role. LGO has a 
pseudocubic structure with an experimental lattice constant of 3.874 A, whereas STO has cubic structure with an 
experimental lattice constant of 3.905 A. We model the Ti0 2 /LaO interface by means of a supercell approach, 
using the average value of these two lattice constants, i.e., 3.884 A. We start from the 6 unit cells of STO, 
containing 120 atoms, and stack 4 unit cells of LGO, contain 80 atoms, on top along the (001) direction. After 
the two slabs we place along the (001) direction a vacuum slab of 15 A thickness. This gives rise to an asymmetric 
model. The supercell contains a total of 200 atoms, of which 16, 24, 16, 24, and 120 are La, Sr, Ga, Ti, and O, 
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respectively. A Ti0 2 layer with four Ti atoms forms the interface to 
the first LaO layer of the LGO slab. Besides the asymmetric model, we 
also address a symmetric model with a total of 292 atoms (32 La, 24 
Sr, 32 Ga, 28 Ti, and 176 O), for which we again apply a vacuum slab 
of 15 A thickness. If not stated otherwise, the following results are 
obtained in both the symmetric and asymmetric models. It has been 
found experimentally that the roughness of the interface (i.e, the 
depth of the cation intermixing) is approximately 2 unit cells for 
the n-type LaA10 3 /STO and LaV0 3 /STO interfaces 5 ' 7 - 27 ' 28 , whereas 
according to Ref. 18 for the LaA10 3 /STO interface a deeper inter- 
mixing might affect the LaA10 3 region. 

The energy required for or gained by cation intermixing is calcu- 
lated as A£ = Eideal — ^intermixed) where Eideal i s me energy of the ideal 
interface and -Eintermixed is the total energy obtained for a supercell in 
which cations have been exchanged. A positive value of AE corre- 
sponds to a stable state after intermixing. Schematic representations 
of the ideal interface structure and different configurations with 
cation intermixing are shown in Fig. 1. In each case, 25% of the atoms 
in the layers affected by intermixing are exchanged. In Fig. 1 we show 
only configurations in which the intermixing is restricted to the 
atomic layers directly at the interface, while we will study in the 
following also diffusion into layers further off the interface. Since 
the lattice relaxation is essential at disordered interfaces, a full relaxa- 
tion of all supercells is required. 

It turns out that the stability (energy gain or loss) of a specific 
intermixing depends on the positions of the exchanged impurity 
atoms. For the TioGa intermixing, we first exchange one interface 
Ti atom with a Ga atom in the 1st, 2nd, 3rd, or 4th LGO unit cell 
(counted from the interface) and then exchange one interface Ga 
atom with a Ti atom in the 1st, 2nd, 3rd, or 4th STO unit cell. The 
same exchange scheme is applied for the SroLa and SrTioLaGa 
intermixing. The calculated values of A£ for TioGa intermixing are 
summarized in Fig. 2. Figure 2(a) shows that when an interface Ti 
atom is exchanged with Ga in the 1st LGO unit cell the system is 
stabilized by 0.06 eV. Similarly, the energy gains for exchange with 
Ga in the 2nd, 3rd, and 4th LGO unit cell are 0.21 eV, 0.29 eV, and 
1.24 eV, respectively. The longer the distance to the interface the 
higher is therefore the energy gain. A similar pattern of the energy 
gain has been obtained for the LaA10 3 /STO interface 4 , because the 
intermixing induces an electric field which compensates the interface 
dipole and thus reduces the energy of the system 17 . This effect is 
demonstrated in Fig. 3, in which we compare the electrostatic poten- 
tials obtained for the ideal and intermixed interfaces. According to 
Fig. 2(b), exchange of an interface Ga atom with Ti in the 2nd, 3rd, 
and 4th STO unit cell destabilizes the system. The energy losses 
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Figure 1 | Structural configurations considered at the LGO/STO 
interface, from left to right: ideal, TioGa intermixing, SroLa 
intermixing, and SrTioLaGa intermixing, in each case directly at the 
interface. 




Figure 2 | Total energy difference induced by the cation intermixing: 
(a) Exchange of an interface Ti atom with a Ga atom in the nth LGO unit 
cell, (b) exchange of an interface Ga atom with a Ti atom in the u-th STO 
unit cell, (c) exchange of an interface Sr atom with a La atom in the u-th 
LGO unit cell, (d) exchange of an interface La atom with a Sr atom in the 
fl-th STO unit cell, (e) exchange of an interface SrTi pair with a LaGa pair 
in the n th LGO unit cell, and (f) exchange of an interface LaGa pair with a 
SrTi pair in the nth STO unit cell. By n = 0 we denote the ideal unit cell. 

amount to -5.99 eV, -4.25 eV, and -3.12 eV. While Ti diffusion 
into the LGO thus is energetically favorable, Ga diffusion into the 
STO is not. The authors of Ref. 22 report a Ti signal past the interface 
on the LGO side and note that this could be due to Ti diffusion. While 
this scenario would agree with our findings, those authors believe 
that the signal is rather due to pollution introduced during the sam- 
ple preparation. 

Figure 2(c) deals with SroLa intermixing. Exchange of an inter- 
face Sr atom with La in the 1st LGO unit cell leads to an energy gain of 
0.05 eV. A similar stability is found for La/Sr disorder at the n-type 
LaTi0 3 /STO and LaA10 3 /STO interfaces 19,29 . For deeper diffusion of 
Sr into the LGO, the total energy is higher than that of the ideal 
system. For A£ we obtain —1.45 eV, —2.32 eV, and —2.83 eV for 
incorporation in the 2nd, 3rd, and 4th LGO unit cell, respectively. 
Moreover, exchange of an interface La atom with a Sr atom in the 
2nd, 3rd, and 4th STO unit cell leads to energy losses of —3.20 eV, 
— 2.59 eV, and —2.48 eV, respectively. We obtain the most stable 
state when La and Sr are exchanged directly at the interface. While La 
diffuses further into the STO, Sr diffusion into the LGO is not pos- 
sible. The coupled SrTioLaGa cation intermixing, i.e., exchange of 
SrTi and LaGa pairs, is addressed in Figs. 2(e) and 2(f). We obtain an 
energy gain only when the exchange occurs directly at the interface. 
While cooperative diffusion of SrTi into the LGO is possible, dif- 
fusion of LaGa into the STO is energetically not favorable. 

Orbitally resolved Ti 3d densities of states (DOSs) for different 
interface configurations are shown in Fig. 4 to address the influence 
of the cation intermixing on the electronic structure. The top row 
gives results for the ideal interface and for TioGa intermixing in the 
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Asymmetric model 




Symmetric model 
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Figure 3 | Electrostatic potential, refering to the most stable configurations of the (top left) ideal interface, (top right) Ti<s>Ga exchange, (bottom left) 
SroLa exchange, and (bottom right) SrTioLaGa exchange. 



4th LGO unit cell (most stable configuration), respectively. Clearly, 
TioGa intermixing has a significant effect on the interface conduc- 
tivity. The Ti states shift to higher energy and become completely 
depopulated, which gives rise to a large band gap. DOSs obtained for 
SroLa and SrTioLaGa intermixing (most stable configurations) are 
shown in the bottom row of Fig. 4. They reveal almost the same 
features as obtained for TioGa intermixing, but with smaller band 
gaps. 

In general, the crystal field due to the octahedral coordination by O 
atoms splits the Ti 3d states into high energy e g (d 3z 2_ r 2 and d x i_ y i) 
and low energy t 2g (d xy , d xz , and d yz ) states. Note that the O octahed- 
ron is distorted at the interface and the crystal field therefore is not 
perfectly octahedral. In case of the clean interface the d xy orbital 
(which is oriented parallel to the interface) is occupied and therefore 
carries the two-dimensional electron gas. The remaining t 2g orbitals 



(d xz and d yz ) are degenerate and occupy the energy range from 
slightly above E P up to about 3 eV. The e g orbitals stay far above 
E F and therefore play no role for the interface metallicity. The charge 
carrier density in the d xy orbitals amounts to 3.4 ■ 10 13 cm -2 , 1.6 ■ 
10" cm" 2 , and 0.1 ■ 10 13 cm" 2 in the 1st, 2nd, and 3rd Ti0 2 layer 
from the interface. Rather similar shapes of the orbitally resolved 
DOSs are obtained for TioGa intermixing, for SroLa intermixing, 
and for SrTioLaGa intermixing, see Fig. 4. However, in each of these 
cases the states appear at much higher energy and no metallicity is 
induced. 



Discussion 

Let us now turn to the question why cation intermixing in LaGa0 3 / 
STO results in a suppression of the interface metallicity, while in 
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Figure 4 | Orbitally resolved partial 3d DOS of the interface Ti atom next to the site where an atom has been exchanged, refering to the most stable 
configurations of the (top left) ideal interface, (top right) TioGa exchange, (bottom left) SroLa exchange, and (bottom right) SrTioLaGa exchange. 



LaTi0 3 /STO no such effect is found 19 . We focus on the first Ti0 2 
layer directly at the interface. As Ti-O-Ti bonding is replaced by 
Ti-O-Ga bonding, differences in the electronegativity between Ti 
and Ga are expected to have important implications. A Ga atom 
contributes less charge to the covalent Ga-O bond due to its high 
electronegativity of 1.81, as compared to the Ti electronegativity of 
1.54. As a consequence, the orbital overlap between Ti and O is 
reduced and delocalized states within the modified Ti0 2 layer 
become energetically less favorable (the band width of the d xy states 
decreases). In addition, the ionic nature of the bonds will also influ- 
ence neighbouring atoms so that a rather small amount of impurities 
can strongly counteract the creation of the two-dimensional electron 
gas. In the case of SroLa intermixing we have to take into account 
that the electronegativity of Sr (0.95) is smaller than that of La (1.1), 



although the difference of 0.15 between Sr and La is smaller than the 
difference of 0.27 between Ti and Ga. This effect can clearly be seen in 
the DOS. The band gap for TioGa intermixing is larger than for 
SroLa intermixing, which confirms that the cation electronegativity 
plays an important role for the two-dimensional electron gas at 
pervoskite oxide interfaces. Therefore, the suppression of the two- 
dimensional electron gas by cation intermixing is a consequence of 
an enhanced ionic character of the metal-0 bonds and the induced 
electric field that compensates the interface dipole. 

In conclusion, we have studied the stability and electronic struc- 
ture of the rc-type LGO/STO interface by first principles calculations. 
We find a two-dimensional electron gas for the ideal interface, while 
the interfaces for TioGa, SroLa, and SrTioLaGa intermixing 
exhibit insulating states. All three types of intermixing can result in 
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an energy gain and therefore in thermodynamic stability. The energy 
gain or loss due to the cation intermixing depends strongly on the 
distance of the exchanged ions from the interface. Importantly, the 
interface electronic structures of the cation intermixed configura- 
tions are all characterized by a complete suppression of the two- 
dimensional electron gas present at the clean interface. This fact 
demonstrates the key role of an atomically clean interface to obtain 
a two-dimensional electron gas. 

Methods 

Our calculations are performed in the framework of spin- degenerate density func- 
tional theory using the generalized gradient approximation in the Perdew-Wang 
flavor and projector augmented wave pseudopotentials, as implemented in the 
Vienna Ab-initio Simulation Package. The presented results have been obtained 
without onsite Coulomb interaction. However, we have checked the metallic cases for 
an onsite interaction of U = 4 eV on the Ti 3d orbitals and find no qualitative 
difference. Relativistic effects are taken into account fully for the core states, while the 
scalar relativistic approximation is used for the valence states (i.e., spin-orbit coupling 
is neglected). The electronic wave function is expanded with a kinetic energy cutoff of 
400 eV. We optimize the crystal structures by minimizing the atomic forces until all 
residual forces remain below 0.001 eV/A. A 10 X 10 X 1 fc-space grid, comprising 21 
points in the irreducible wedge of the Brillouin zone, is found to be well converged. 
Self-consistency is assumed for a total energy convergence below 0.0001 eV. 
Moreover, the DOS is calculated with a Gaussian smearing of 0.05 eV. 
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